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ABSTRACT 

Dust grains coagulate into larger aggregates in dense gas. This changes their size distribution 
and possibly affects the thermal evolution of star-forming clouds. We here investigate dust co- 
agulation in collapsing pre-stellar cores with different metallicities by considering the thermal 
motions of grains. We show that coagulation does occur even at low metallicity ~ 10~ 6 Zq. 
However, we also find (i) that the H 2 formation rate on dust grains is reduced only after the 
majority of H2 is formed; and (ii) that the dust opacity is modified only after the core becomes 
optically thick. Therefore, we conclude that the effects of dust coagulation can safely be ne- 
glected in discussing the temperature evolution of the pre-stellar cores for any metallicity as 
long as the grain motions are thermal. 

Key words: dust, extinction — galaxies: evolution — galaxies: high-redshift — ISM: abun- 
dances — molecular processes — stars: formation 



1 INTRODUCTION 

The first gene ration of stars is believe d to be very massive (> 
100 M ; e.g.i romm & La rson 2004), while those in the sola r 
neighbourhood have typical masses of < M (e.g. lLarsonll2005h . 
The transition in the characteristic stellar mass should have oc- 
curred in the history of the Universe. The most widely accepted 
cause fo r this transition is metal enrichment in the interstella r 
medium dOmukail200b1:lBromm et alj200l|jschneider et alj2002h . 
In other words, the transition of the characteristic stellar mass oc- 
curs at a 'critical metallicity'. 

Dust grains should be produced in s upernovae of first stars 
and thus be present in metal-enriched gas ( Todini & Ferrara 2001; 
iNozawa et atll2003l : ISchneider. Ferrara. & Salvaterrall2004l) ! Their 
thermal emission becomes effective for gas cooling through colli- 
sional coupling at such high density that the leans mass is below 
Mq. With a sufficient amount of dust, dust cooling is able to in- 
duce fragmentation into subsolar mass cores. If most of the metals 
condense into dust, t he critical metallicity is as low as ~ 10 Zq 
dOmukai et"ai]2 005), although difference in grain composition and 
grain size d istribution alters the exact value by about an order of 
magnitude dSchneider e t al. 2006). 

In such dense clumps where dust cooling is efficient, how- 
ever, coagulation of dust grains could proceed and modify the dust 
opacity as well as the H2 formation efficiency on dust grains, con- 
sequently affecting the thermal evolution of star-forming clouds. 
Therefore, coagulation may have a large influence on the critical 
metallicity. In spite of this potential importance, the effects of co- 
agulation on the thermal evolution of star-forming clouds have not 
been considered in previous studies. 
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In this paper, we assess the effects of coagulation on the ther- 
mal evolution of star-forming clouds. We focus on low metallicity 
in the context of the critical metallicity, but we treat all the range 
of metallicity from nearly zero to Zq . Therefore, the methods and 
results in this paper are applicable to star-forming clouds at any 
metallicity. This paper is organized as follows. The method used in 
this paper is explained in Section[2] We describe some basic prop- 
erties of the thermal evolution of collapsing cores in Section[3] The 
effects of coagulation are discussed in Section[4] Finally, Section|5] 
gives the conclusion. 



2 METHOD 

We examine the effect of coagulation on the thermal evolution of 
star-forming collapsing cores. First, the density and temperature 
evolution of collapsing cores are calculated with a fixed grain size 
distribution (i.e. without coagulation). Then, under these calculated 
density and temperature, we estimate the effect of dust coagulation. 
We explain the basic ingredients of our calculations in the following 
with particular emphasis on physics related to dust grains. 



2.1 Basic equations for pre-stellar collapse 

We calculate the thermal evoluti on of star-form ing cores by using 
the formulations dev eloped by lOmukail d2000h and improved by 
lOmukai et al.l (12005). Here we briefly review the methods. 

The collapsing gas in a star-forming core is treated as one zone 
with a homogeneous hydrogen nuclei number density nu and a 
single gas temperature T. These density and temperature can be 
interpreted as those in the central region of the core, whose size is 
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given by the Jeans length. The evolution of density is calculated by 

P .., I 3-7T 



dp 
dt 



with tff = 



32Gp' 



(1) 



where ts is the free-fall time and G is the gravitational constant. 
The temperature evolution is calculated by solving the energy equa- 
tion, 



de 



dt P dt \ p 



An 



(2) 



where e is the specific thermal energy, p is the pressure, p is the 
density, and A ne t is the net cooling rate. The net cooling rate in- 
cludes all the important processes: molecular/atomic line emission, 
dust and gas continuum emission, and chemical heating and cool- 
ing. For the details o f the coo ling/heating processes and of their 
treatment, we refer to Om ukail feOOOl) . although we review the pro- 
cesses related to dust grains in the following subsections. We addi- 
tionally impose a cosmic background radiation at 1 + z = 10 (i.e. 
Tomb = 27 K), so that the temperature does not drop below it. 
However, our results on coagulation are not sensitive to the back- 
ground temperature. 



2.2 Dust cooling 

Dust grains achieve the following balance between dust thermal 
emission and heating due to collisions with gas particles: 



4a(T gr - TcmbKp /?escpl> = H gI , 



(3) 



where a is the Stefan-Boltzmann constant, T gI is the dust temper- 
ature, ftp is the Planck mean opacity of dust grains per unit dust 
mass, /3 CS c is the photon escape probability, T> is the dust-to-gas 
mass ratio, and H gT is the collisional he ating (cooling) rate of dust 
(gas) particles. Following Omukai (2000), the photon escape prob- 
ability is written as f3 CBC — min(l, t -2 ) by using the optical depth 

T. 

The collisional heating rate of dust particles can be written as 
dHollenbach & McKeell 19791) 



Ha 



n gr (2fc B T- 2fc B T gr ) 
icon 



(4) 



where t~ oll = nna gT VHf is the mean interval between two suc- 
cessive collisions, n gI and <r gr are the grain number density and 
cross-section, respectively, vh = (8ksT /tyitiu) 1 ^ 2 is the aver- 
age speed of hydrogen nuclei under a Maxwellian distribution, and 
/ = 0.3536 + 0.5yHc (j/He = 0.083, corresponding to a He mass 
fraction of Y = 0.25) measures the contribution from species other 
than protons. Once the grain temperature exceeds the sublimation 
temperature, T su b, the grains are assumed to disappear. 

The grain cross-section is related to the total surface area of 
grains per unit dust mass, S, as 



rmmH(l + iync)ST>. 



(5) 



Given the grain size distribution functions per unit dust mass, 
-<V g r,i(a) (i.e. Af S i,i(a) da gives the number of dust grains of 
species i with radius between a and a + da per unit dust mass 
and the grains are assumed to be spherical), and the mass fractional 
abundance, fi, of each grain species i, S is quantified as 



S = ^2fiSi with S, 



Af g r,i(a)na 2 da. 



(6) 



surface, and the line emission from formed H2 contributes to gas 
cooling. 

The total absorption coefficient at frequency v per unit dust 
mass is written in a similar way as 



k{v) = Y^/i/t^i with 



>-;, , = I Q % v (a)N g ? ,i(a)na da, (7) 




where Q„(a) is the absorption cross-section of species i normal- 
ized to the geometrical cross-section. 



2.3 H2 formation on dust grains 

The H2 formation rate on grain surface per unit volume can be 
expressed as Rn 2 = k gr n(H)nu, where n(H) is the number 
density of H atoms in the g as phase. We use the rate coefficient 
dTielens & Hollenbach|[l985 



fc gr = 6 x lO- 17 (S'/S ) v /r/3OOK/ a 5*H, 



(8) 



The grain surface area determines the H2 formation rate on grain 



where S is the total surface area of dust grains, Sq is that for the 
MRN grain size distribution (Section [2.51 ) for Z = Zq (Z is the 
metallicity), 

/a = {l + exp [7.5 x 10 2 (l/75 - 1/T gr )] } _1 , (9) 
and 

Sh = [1 + 0.04(T + T gr ) ' 5 + 0.002T + 8 x 10~ 6 T 2 ] -1 . (10) 

Note that, without coagulation, S/Sq is simply Z/Zq. The rate 
coefficient (equation [8) with S/Sq — 1 is consistent with obser- 
vationally derived rates for a sample of clouds in the solar neigh- 
bourhood jjurall975l : lHollenbach & McKedl979h . Thus we expect 
that equation {8]l with S/ Sq — 1 is valid for the MRN grain size 
distribution and solar metallicity. 



2.4 Grain coagulation 

The time evolution of the grain size distrib ution by coagulation i s 
calculated by adopting the formulation of iHirashita & Yar] J2009h . 
Since they treated coagulation under a constant density, we mod- 
ified the formulation to treat the varying gas density in collapsing 
cores. We briefly review the calcul ation method. The details are de- 
scribed in Hiras hita & Yar] d2009h . In this subsection, we omit the 
index i specifying t he species. 

As assumed in IHirashita & Yar] ( |2009h . we consider grains to 
be compact and spherical. In other words, we do not consider fluffy 
aggre gates, which are suggested to form after sticking of grains 
(Meak in & DonnlligSSlJOssenkoDjliggBh Fl There are two compet- 
ing effects of the fluffi ness. Fluffy grains wo uld increase the coagu- 
lation rate l lOrmel. Spaans, & Tielens 2007), which could enhance 
the decrease of S. On the other hand, fluffy grains have larger S 
than compact grains. Because of these competing effects on S, it 
is not obvious if the change of S is under- or overestimated by as- 
suming a compact structure. 

We solve the coagulation equation discretized for the grain 
size. The total grain mass density is normalized to nnmn(l + 



We note that aggregation rather than coagulation may be more appro- 
priate terminology to express the grain growth after grain-grain collisions, 
especially when we consider fluffy aggregates. However, we use the term 
'coagulation' throughout this paper, since we only consider compact (i.e. 
not fluffy) grains. 
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4yHc)2?- Thus, as tih becomes denser in the evolution of collapsing 
cores, the grain density also becomes proportionally larger. 

We consider silicate and graphite as grain species. In order 
to avoid complexity in compound species, we only treat collisions 
between the same species. The time evolution of p k (grain density 
contained in the fc-th size bin; all the subscripts in the following 
equation specify the bin number) by coagulation can be written as 



dpfc _ 
dt 

where 



-m k p k 



1\ iV iV 

^a ek p e + -^2^2a ej p e pjm e c 3 oag (k), (11) 



J if «tt < Wcoag> 

O-lk = < m k me 

otherwise, 



(12) 



(f) is the sticking probability) and m l J os _ z {k) = nit if the mass 
of the coagulated particle mi + rrij is in the mass range of the 
fc-th bin; otherwise m? oag (l:) = 0. Coagulation is assumed to 
occur only if the relative velocity is less than the coagulation 
threshold velocity «coag- The cross-section for the coagulation is 
ffifc = Tt(at + a k ) 2 . The coagulation th r eshold velocity is given by 



dChokshi Tielens. & Hollenbacrj[l993l ; iDominik & Tielens|[l99 : f 
Yan, Lazarian, & Draine 2004) 



21.4 



at + a! 
(a e + a k y 



1/2 



v 5/6 



Rl/3 ft ' n 



5/6 1/2' 



(13) 



where 7 is the surface energy per unit area, Ri k = aeak/(a-e + a fe) 
is the reduced radius of the grains, E is related to Poisson's ra- 
tios (yi and u k ) and Young's modulus (Et and E k ) by 1/E = 
(1 - Vtf/E t + (1 - v k f/E k . We also assume (3 = 1 for the 
sticking probability. A sticking probability of order unity below 
the co agulation threshold is shown by experimental studies (Blum 
2000). The effect of /3 is also addressed in Section PTTl 

We adopt a single velocity for each grain size. The velocity is 
given by the thermal (Brownian) velocity as 



/ 8fc B T \ 
V nm k J 



1/2 



(14) 



Smaller grains have larger thermal velocities, and the effect of co- 
agulation firs t appears fo r the sm allest grains (see also Section[4]l. 
As shown in Ossenkopl dl993h . motions other than the thermal 
motion are potentially important. In particular, large grains tend 
to be coupled with turbulent motion if it exists ( Ossenkopfl l 19931 ; 
IWeidenschilling & Ruzmaikinall 19941 ; lOrmel et alj|2009l) l 2 l For the 
sake of simplicity and exploratory nature of this work, the contri- 
bution of other velocity sources is ignored here and remains open 
for further study. 

Each time-step is divided into four equal small steps, and 
we apply v k £ = v k + ve , \v k — Ve\, v k , and v% in each step, 
following Hira shita & Yar] d2009l) . A relative velocity of v k t = 
(SksT/np) 1 ^ 2 , where p the reduced mass, is often adopted, but 
the difference caused by the different treatment of relative velocity 
is too small to change our conclusion. 



2 Recently. lOrmel et alj J20091) show that even if turbulent velocity is con- 
sidered, large grains with a > 0.1 pm do not coagulate significantly on 
the free-fall timescale. However, they consider relatively small densities 
(< 10 7 cm -3 ), and the effect of turbulence on the higher density range 
is still open. 



2.5 Adopted parameters for dust grains 

As stated above, we assume a mixture of silicate and graphite 
species. We adopt T> — 0.006(Z/Zq) for both silicate and 
graphite, and the contribution from these two species are summed 
with weights ft = 0.54 and 0.46 for silicate and graphite, respec- 
tively dHirashita & Yarl2009l) . The dust-t o-gas rat i o in so lar metal- 
licity environments (0.006) is take n from Spitzei dl978 ). The ma- 
terial densities s are adopted from lDraine & Led dl984l) (3.3 and 
2.26 g cm -3 for silicate and graphite, respectively). 

The quantities for c o agulati on are taken from 
Choks hi. Tielens. & Hollenbachj [l993h (the data for quartz 
and graphite are used for silicate and graphite, respectively; i.e. 
7 = 25 and 12 erg cm" 2 , E = 5.4 x 10 11 and 3.4 x 10 10 
dyn cm -2 , and v — 0.17 and 0.5 for silicate and graphite, 
respectively). The sublimation temperatures (T su b = 1400 and 
1750 K for si l icate and graphite, respectively) are taken from 
lLaor& Draine! dl993h . 

The absorption cross-section of grains is calculated by us- 
ing the Mie theory dBohren & Huffmanlll983h . The optical con- 
stants of silica t e (astr onomical silicate) and graphite are taken from 
iDraine & Leel {j984). The grain size distribution is assumed to be 
A/" gr ,i(a) oc a~ 3 ' 5 in the size range 0.005 pm ^ a ^ 0.25 pm 
and A/" gr .i(a ) = out of this range for both silica te and graphite as 
indicated by Mat his. Rumpl. & Nordsieckl dl977h (MRN). 



3 THERMAL EVOLUTION WITHOUT COAGULATION 

First, density and temperature evolutions are calculated without co- 
agulation (i.e. with the fixed MRN grain size distribution; Section 
12.5b - Then, the density and temperature given in this section are 
used as the background density and temperature under which grain 
coagulation is calculated in the next section. 

The evolutions of collapsing cores are shown in Fig. Q] in the 
metallicity range 10" 7 < Z/Zq 1. In the case with 10~ 7 Zq, 
metal effects are negligible and its temperature evolution is identi- 
cal to the metal-free one. In cores with metallicity < 10~ 4 Zq, H2 
is mainly formed via gas-phase reactions. With higher metallicity, 
i.e. higher dust abundance, the majority of H2 is formed via the dust 
surface reactions. H2 line cooling contributes significantly to the 
thermal evolution at relatively low densities («h < 10 5-8 cm -3 ) as 
seen by the dips of the temperature evolution in this density range, 
although carbon and oxygen fine-structure lines contribute more 
to the cooling at Z > 10~ 2 Zq. Therefore, the modification of 
grain surface area (if any) is the most relevant around this density 
range. On the other hand, the impact of dust cooling appears around 
n H ~ 10 10 " 13 cm -3 for metallicities Z > 10" 6 Zq as seen by the 
dips of the temperature evolution. At higher densities, dust cooling 
has minor influence on the thermal evolution, since the system be- 
comes optically thick to dust continuum. Although the metallicity 
range wh ere the effect o f dust cooling appears is broadly consis- 
tent with lOmukai et at] d2005h . the effect is more pronounced at 
the same metallicity in our case, especially for the lowest values 
~ 10 -6 Zq. This is because our assumed composition of carbona- 
ceous grai ns, graphite, is more r efractory than organic compounds 
adopted in Omu kai et al.1 d2005t) . Grains survive up to nn ~ 10 15 
cm -3 in our case. Our choice is conservative in the sense that the 
coagulation effect becomes more visible as the dust survives and 
coagulates up to higher temperature and density. 
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Figure 1. Temperature evolution of collapsing pre-stellar cores with dif- 
ferent metallicities as a function of the hydrogen number density, which 
increases with time. Solid curves indicate the cases with metallicities of 
log(Z/Z0) = —7, —5, —3, and —1, while dashed curves show the re- 
sults with log(Z/Zo) = —6, —4, —2, and 0. The numbers represent the 
values of log(Z/Z0) adopted for individual lines. 

4 EFFECTS OF COAGULATION 

Here we evaluate the effect of coagulation under the density- 
temperature evolution calculated in the previous section. Coagula- 
tion could affect the thermal evolution of collapsing cores through 
the changes of S and ftp. The former and latter quantities have 
influence on the H2 formation rate and the dust cooling rate, re- 
spectively. In Fig. [2] we show the evolution of these two quan- 
tities. Both quantities are divided by the values without coagu- 
lation (i.e. the MRN grain size distribution) and are denoted as 
K P (coag)//tp(MRN) and S'(coag)/S(MRN). The effects of co- 
agulation appear if these ratios deviate significantly from unity. In 
fact, they change at high density, which indicates that coagulation 
really takes place efficiently enough to affect the grain opacity and 
the grain surface area. 

In Fig. [2] we also show the dust cooling rate relative to 
the total (compressional and chemical) heating rate and the in- 
crease rate of H2 fraction per free-fall time. The latter is equal to 
(d/H 2 /dlnnH)dust, where /h 2 is the H2 fraction (/h 2 = 1 if hy- 
drogen is fully molecular), and the subscript "dust" indicates the 
contribution from grain surface reaction. 

In order to have an idea about how coagulation proceeds, we 
also show the grain size distributions in Fig. [3] where the grain size 
distribution is multiplied by a 4 to show the grain mass distribution 
per logarithmic size (the largest contribution to the dust mass comes 
from the grain size range where a 4 J\f gT peaks). By comparing Fig. 
[3] with Fig. [2] we observe that the decrease of grains with a ~ 
ctmin is indeed seen at the density range where 5(coag)/5'(MRN) 
decreases in Fig. [2] The increase of Kp(coag)/Kp(MRN), on the 
other hand, occurs when the grains are coagulated to a > 10~ 5 cm. 
In the following, we will discuss further details. 

4.1 Grain surface area and coagulation criterion 

We compare the dashed lines of the upper and lower panels (i.e. 
grain surface area vs. H2 formation rate by grain surface reaction) 
in Fig. [2] We observe that the drop of grain surface area generally 
occurs well after the H2 formation on grain surface is completed. 
This means that coagulation has no influence on the thermal evo- 



lution of collapsing cores either through H2 line cooling or H2 for- 
mation heating. The main reason why H2 formation is faster than 
grain coagulation is that the thermal velocity of H atoms is much 
larger than that of grains. 

The coagulation time-scale t CO ag can be estimated by the 
grain-grain collision time-scale. The condition that coagulation 
takes place on a free-fall time-scale (t CO ag < iff) is expressed as 
n-R > 7iH,coa g , where jiH, C oag is called "critical density for coag- 
ulation" in this paper. The critical density for coagulation is esti- 
mated as 



^H,coag 




where s is the grain material density. This equation is derived in 
the Appendix, where the grain size is represented by a single value 
(a) for the simplicity of estimation. Indeed as shown in Figs.|2]and 
[3] if tih > iiH.coag, coagulation significantly decreases the grain 
surface area, and the grain size distributions at a < 10~ 6 cm are 
indeed affected by coagulation. This criterion can be used gener- 
ally, as long as a grain size distribution similar to MRN is adopted 
and the grain velocity is thermal. Equation (JT3J can also be solved 
for a under a given hydrogen number density to obtain a condition 
for the grain size affected by coagulation. 

As shown in equation J 15b . if we assume a smaller sticking ef- 
ficiency (J3 < 1), riH.coag increases in proportion to [3~ 2 . Thus, the 
conclusion that coagulation has no influence on the thermal evolu- 
tion of collapsing cores is strengthened if we adopt < 1. 

4.2 Grain opacity 

We compare the solid lines of the upper and lower panels (i.e. 
Planck mean opacity vs. dust cooling) in Fig. [2] We find that 
the change of Kp (coag) /kp (MRN) appears after the epoch in 
which dust cooling dominates. In other words, before coagula- 
tion affects the grain opacity, the contribution from dust cool- 
ing has already declined because the system becomes optically 
thick to dust continuum. This indicates that the change of grain 
opacity by coagulation has no influence on the thermal evolution. 
In particular, Kp (coag) / Kp (MRN) is kept to be unity even after 
S'(coag) /S'(MRN) drops significantly. This is because Kp does not 
depend on the grain radius under a fixed total grain mass. In fact, 
from equation 10, Kp oc a Qi/Ngr (A^g r is the total grain number 
per unit dust mass) with a e valuated at a typical g rain radius, while 
Q v oc a if A = c/v > a jDraine & Lea l "l 984) and N gI oc a~ 3 
from the conservation of the total dust mass. 

Since the grain opacity is proportional to the grain mass, Fig.|3] 
is useful to understand which grain size dominates the grain opac- 
ity. The largest grains dominate the grain mass in the MRN grain 
size distribution, which means that the grain opacity is also domi- 
nated by the largest grains. 

In all the metallicities except for Z = 10~ 6 Zq, we find 
that Kp (coag) / Kp (MRN) increases at a high density depending on 
metallicity. At this density, the dust temperature, which is almost 
equal to the gas temperature, is ~ 100-1000 K, and the radiation 
spectrum has a peak around wavelengths of A ~ 3-30 /im. At this 
wavelength range, Q v /a of graphite increases as a function of a if 
a > 0. 1 /xm because of the contribution from magnetic dipole ra- 
diation jDraine & Leell 19841) . Moreover, the fraction of grains with 
a > 0.1 fim is enhanced after coagulation. These two effects (tem- 
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Figure 2. Upper panels: the evolutions of Planck mean dust opacity (solid line) and grain surface area (dashed line). The values of these two quantities in the 
presence of coagulation are normalized to those without coagulation (i.e. with the initial MRN grain size distribution preserved). Lower panels: the ratio of 
dust cooling to the total heating rate (solid line), and the increase rate of H2 fraction by the grain surface reaction in a free-fall time (dashed line). Each panel 
shows the result with metallicity (a) 10~ 6 Zq, (b) 1CT 4 Zq, (c) 10" 2 Z , and (d) 1 Z Q . 



perature increase and grain growth) are the reason for the increase 
of k p (coag) /k p (MRN). 



5 CONCLUSION 

By calculating the thermal evolution and the grain coagulation in 
collapsing cores, we have shown that coagulation does not affect 
their thermal evolution for any metallicity (0 ^ Z < 1 Zq) as 
long as the grain velocities are governed by thermal motions. 

Although the above statement is the most important conclu- 
sion in this paper, there are two things that are worth raising. First, 
coagulation really takes place when the density is high enough. 
Coagulation is important at rtn > 10 7 cm -3 for Z — Zq and 



the density range with efficient coagulation shifts to higher den- 
sities because of lower grain abundance (i.e. lower grain colli- 
sion rate). The condition that iff > t CO ag is equivalent to hh > 
10 7 (Z/Zq)~ 2 (T/100 K)" 1 cm" 3 . This can be widely used as 
the criterion that coagulation occurs in a free-fall time as long as 
the grain size distribution is MRN and the grain velocities are dom- 
inated by the thermal motions. Second, Kp is unaffected by coagu- 
lation up to high densities where the contribution from dust cooling 
drops. The reason why Kp is insensitive to coagulation is that Kp 
is almost independent of the grain radius if a <C A (A is the wave- 
length at which the radiation spectrum peaks) holds and the dust 
temperature is < 100 K. 
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Figure 3. Evolutions of grain size distribution as a function of hydrogen number density, which traces the stage of gas collapse. The grain size distribution 
function per unit dust mass, A/gr(a), is multiplied by a 4 to show the grain mass distribution function. Panels (a), (b), (c), and (d) show the results with 
Z = 10~ 6 , 10~ 4 , 10~ 2 , and 1 Zq, respectively. In each panel, the dotted line presents the initial MRN size distribution, while the solid, dashed, and 
dot-dashed lines represent various evolutionary stages at the densities indicated in each panel. 



ACKNOWLEDGMENTS 

We thank the anonymous referee for useful and constructive com- 
ments which improved this paper considerably. We thank A. Fer- 
rara for helpful suggestions. This study is supported in part by the 
Grants-in-Aid by the Ministry of Education, Culture, and Science 
of Japan (19047004, 21684007; KO). 



APPENDIX A: DERIVATION OF EQUATION QD 

We provide an estimate of coagulation time-scale. The coagulation 
time-scale, t CO a g , can be estimated by the time-scale on which a 
grain encounter another grain: 

icoag — -z T~Td ' (Al) 

pn gI Tv(a)' : v gI 

where we have assumed that the grains have typical representative 
radius (a.) and velocity « gr , and n gr is the grain number density. 
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We explicitly write the sticking probability /3, which is assumed to 
be 1 in the text. The grain velocity can be estimated by the thermal 
velocity, 



1/2 



7rm gr 

where the typical grain mass m gr is estimated to be 

4tt, , 3 

TO gr ~ —(a) S. 



(A2) 



(A3) 



The grain number density, n gr can be related to the dust-to-gas ratio 
as 



4 o 

-7r(a) sn gr = n H m H (l + 4y Ho )23. 



(A4) 



Using equations l lA2t . iA3l , and l lA4t . equation JAU is written as 

2V67T <a> 5 /V/ 2 



9 /3n H m H (l + 4j/Hc)©(fc B r) 1 /2 



5.7 x 10 1 



"ii 



-/T 1 



(a) 



10" 



23 



,6 x 10" 3 



T 



100 K 



5/2 



-1/2 



3 g cm" 



3/2 



(A5) 



Now we examine if coagulation takes place within the collaps- 
ing timescale, i.e. the free-fall timescale. The free-fall time-scale is 
numerically estimated as 



tff : 



1.4 x 10 1 



(A6) 



Significant coagulation occurs if t C oa g < tg is satisfied, which is 
equivalent to the condition nu > iiH.coag. The critical density for 
coagulation, iin.coag, is given in equation J 15l >. where we have re- 
placed 23/6 x 10~ 3 with Z/Zq by assuming that dust-to-gas ratio 
is proportional to metallicity. 
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